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ABSTRACT 



Aims. This work aims at investigating the molecular gas component in the vicinity of 

two young stellar object (YSO) candidates identified at the border of the HII region 

G034.8-0.7 that is evolving within a molecular cloud shocked by the SNR W44. The 

CO ■ purpose is to explore signatures of star forming activity in this complex region. 

CO ■ 

. Methods. We performed a near and mid infrared study towards the border of the HII 

\ region G034.8-0.7 and observed a 90"x 90" region near 18 h 56 m 48 s , +01°18'45"(J2000) 

' using the Atacama Submillimeter Telescope Experiment (ASTE) in the 12 CO J=3-2, 

13 CO J=3-2, HCO+ J=4-3 and CS J=7-6 lines with an angular resolution of 22". 
Results. Based on the infrared study we propose that the source 2MASS 
18564827+0118471 (IR1 in this work) is a YSO candidate. We discovered a bipolar 
12 CO outflow in the direction of the line of sight and a HCO + clump towards IR1, 
confirming that it is a YSO. From the detection of the CS J=7-6 line we infer the 
presence of high density (>10 7 cm -3 ) and warm (>60 K) gas towards IR1, probably 
belonging to the protostellar envelope where the YSO is forming. We investigated the 
possible genetic connection of IR1 with the SNR and the HII region. By comparing 
the dynamical time of the outflows and the age of the SNR W44, we conclude that 
the possibility of the SNR has triggered the formation of IR1 is unlikely. On the other 
hand, we suggest that the expansion of the HII region G034.8-0.7 is responsible for the 
formation of IR1 through the "collect and collapse" process. 

Key words. ISM: molecules - ISM: clouds - ISM: jests and outflows - stars: formation 
1. Introduction 

The formation of stars occurs deep inside molecular clouds. In the standard scenario for 
the formation of an isolated low-mass star, a cold core contracts as magnetic and turbulent 
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support are lost and subsequently collapses to form a protostar with a surrounding disk. 
Then, a stellar wind breaks out along the rotational axis of the system and drives a bipolar 
outflow, which gradually disperses the protostellar envelope revealing a pre-main sequence 
star with a disk (van Dishoeck & J0rgensen 2008). 

Several authors have studied the role of the supersonic turbulence in the star forming 
processes (e.g. Mac Low & Klessen 2004). They suggest that when stars form the process is 
quick and dynamic, with gravitational collapse occurring at a rate controlled by supersonic 
turbulence. The balance between turbulence and gravity provides a natural explanation 
for the widely varying star formation rates seen at both cloud and galactic scales. 

As the largest contribution to interstellar turbulence comes from supernova explosions, 
it is thought that the dominant driving mechanism in star-forming regions of galaxies 
appears to be supernovac, while elsewhere coupling of rotation to the gas through magnetic 
fields or gravity may be important (Mac Low & Klessen 2004). 

Since the early suggestion made by Opik (1953) several numerical studies (e.g. Vanhala 
& Cameron 1998; Mclioli et al. 2006) proposed that supernova remnants (SNRs) can initiate 
the formation of new generations of stars. However, observational evidence of star formation 
triggered by SNRs is still scarce. The bright SNR W44 and the environment where it is 
evolving is a good target to explore this possibility. W44 is one of the few demonstrated 
cases of a SNR-molecular cloud interacting system. Seta et al. (1998) discovered six giant 
molecular clouds that appear to be partially surrounding the remnant. Later, Seta et al. 
(2004) proved that some of these clouds (at vlsr ~ 45 km s^ 1 ) are physically interacting 
with the remnant on its southeastern side. The CO spectra corresponding to these molecular 
structures present wing emission that, as the authors conclude, unambiguously confirm the 
interaction between W44 and the molecular gas. The physical interaction of the blast wave 
with a clumpy interstellar medium is also supported by the detection of bright OH (1720 
MHz) masers at LSR velocities approximately between 43 and 47 km s^ 1 (Hoffman et al. 
2005, and references therein) and shocked H2 emission observed with Spitzer-IRAC (Reach 
et al. 2006). Many infrared sources are embedded in the portion of the giant molecular 
cloud that is being shocked by the SNR. In particular IRAS 18542+0114, an IRAS point 
source that coincides with the border of the HII region G034.8-0.7 (Kuchar & Clark 1997; 
Ortega et al. 2007) appears immersed in a region of shocked molecular gas interior to the 
SNR shell as seen in projection. Figure 1 shows the molecular cloud in the 13 CO ,1=1 
line (as extracted from the GRS 1 ) from the average between 40 and 50 km s -1 , the velocity 
range where the molecular cloud peaks, and also in coincidence with the systemic velocity 
of W44. The radio continuum of the SNR and the HII region at 20 cm as extracted from 
the MAGPIS (Hclfand et al. 2006) arc depicted in thick blue contours. 

The distance to W44 and the shocked molecular cloud was estimated to be ~ 3 kpc 
(Radhakrishnan et al. 1972; Caswell et al. 1975; Seta et al. 2004). Taking into account 
the recombination line velocity of 52.1±3.6 km s _1 derived by Lockman (1989) and the 
HI absorption study performed by Ortega et al. (2007), the same distance as the SNR is 
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suggested for the HII region G034.8-0.7. In what follows we adopt 3 kpc as the distance for 
the SNR W44, the giant molecular cloud GMC G34.8-0.6 and the HII region G034.8-0.7. 

IRAS 18542+0114 (cross in Figure 1) is resolved into several sources in the Two Micron 
All-Sky Point Source Catalog (2MASS, Cutri et al. 2003), and according to color criteria 
(see Section 3.1) some of them are young stellar object (YSO) candidates. 



Fig. 1. 13 CO J=l-0 line averaged between 40 and 50 km s showing the giant molecular cloud 
GMC G34.8-0.6. The thick blue contours outline the radio continuum emission of the SNR W44 
and the HII region G034.8-0.7 at 20 cm. The cross shows the position of IRAS 18542+0114. 

Thus, the scenario is that there are YSO candidates within a shocked molecular cloud 
located in a region peripheric to a SNR and simultaneously on the border of an expanding 
HII region. 

In this work we present a near and mid infrared study and molecular observations in 
the 12 CO and 13 CO J=3-2, HCO+ J=4~3 and CS J=7-6 lines towards IRAS 18542+0114 
carried out with the purpose of exploring the star forming activity in the region and its 
connection with the local SNR and the HII region. 

2. Observations 

The molecular observations were performed on June 25, 2008 with the 10 m Atacama 
Submillimeter Telescope Experiment (ASTE; Ezawa ct al. 2004). We used the CATS345 
GHz band receiver, which is a two-single band SIS receiver remotely tunable in the LO 
frequency range of 324-372 GHz. We simultaneously observed 12 CO J=3-2 at 345.796 GHz 
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and HCO + J=4 3 at 356.734 GHz, mapping a region of 90" x 90" centered at the position 
of IRAS 18542+0114 (RA = 18 h 56 m 47.8 s , dec. = +01°18'45", J2000). The mapping grid 
spacing was 10" and the integration time was 72 sec. per pointing. Additionally we observed 
13 CO J=3-2 at 330.588 GHz and CS J=7-6 at 342.883 GHz towards the centre of the region. 
All the observations were performed in position switching mode. The off position was RA 
= 18 h 55 m 30 s , dec. = +01°39'55" (J2000) that was checked to be free of emission. 

We used the XF digital spectrometer with a bandwidth and spectral resolution set 
to 128 MHz and 125 kHz, respectively. The velocity resolution was 0.11 km s _1 and the 
half-power beamwidth (HPBW) was 22" at 345 GHz. The system temperature varied from 
T sys = 400 to 700 K. The typical rms noise (in units of T m b) ranged between 0.1 and 0.4 
K and the main beam efficiency was fy m b ~ 0.65. 

The spectra were Hanning smoothed to improve the signal-to-noise ratio and only lin- 
ear or/and some third order polinomia were used for baseline fitting. The spectra were 
processed using the XSpec software package developed at the Onsala Space Observatory. 



3. Results 

3.1. Infrared emission 

Figure 2 shows a Spitzer-IRAC three color image of a region about 22' x 22' around IRAS 
18542+0114 as extracted from GLIMPSE (Fazio ct al. 2004; Werner ct al. 2004). The three 
IR bands arc 3.6 /im (in blue), 4.5 /im (in green) and 8 /im (in red). Contours of the radio 
continuum emission at 20 cm of the SNR W44 and the HII region G034.8-0.7 are included. 

The IRAC 4.5 fan band contains both H 2 (v = 0-0, S(9, 10,11)) lines and CO {v = 1-0) 
band heads. As noticed by Cyganowski et al. (2008), all of these lines may be excited by 
shocks, such as those expected when protostcllar outflows crash into the ambient ISM. In 
Figure 2 it can be noticed that IRAS 18542+0114 appears slightly extended in the 4.5 
fim emission (green), therefore suggesting that it may be a YSO outflow candidate. It is 
important to note the difference between the emission from this YSO candidate and from 
the gas shocked by the SNR, both seen in green. The first one is bright and concentrated 
around IRAS 18542+0114, while the second one is seen as diffuse filaments. From Figure 2 
it is also evident the presence of an infrared dark cloud (IRDC 34.776-0554 as catalogued 
in the Glimpse Dark Cloud Catalog). IRAS 18542+0114 is seen in projection over a border 
of this infrared dark cloud (IRDC). The IRDCs have been shown in recent years to be sites 
of the earliest stages of star formation (Rathborne et al. 2006, 2007). 

Figure 3 shows near infrared (NIR) JHK three-color image in a region of about 
10' x 10' that contains IRDC 34.776-0.554 obtained from the 2MASS Survey (Skrutskie 
et al. 2006). The most reddened NIR sources are indicated with numbers. Among these 
sources it is included the bright 2MASS 18564827+0118471, possible counterpart of IRAS 
18542+0114. 
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Fig. 2. Spitzer-IRAG three color image (3.5 /im = blue, 4.5 /im = green and 8 (im = red) of 
a region about 22' x 22' around IRAS 18542+0114. The contours correspond to the continuum 
emission at 20 cm of the SNR W44 and the HII region G034.8-0.7. 



Fig. 3. 2MASS JHK three-color image towards IRDC 34.776-0.554. The most reddened NIR 
sources are indicated with numbers. IRAS 18542+0114 is resolved into several 2MASS sources. It 
can be noticed that the brightest and the most reddened one is 2MASS 18564827+0118471. 



To look for primary tracers of stellar formation activity in the vicinity of the IRDC, we 
used the 2MASS All-Sky Point Source Catalogue in bands J (1.25 /im), H (1.65 /mi) and K 
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(2.17 /im). We performed photometry in the region shown in Figure 3. Figure 4 shows the 
(H-Ks) versus (J-H) color-color (CC) diagram. We selected sources with Ks < 15 and with 
detection in at least two bands. The typical color errors are about 5 %. Based on the above 
criteria we found 1029 sources in the mentioned region. Following Hanson et al. (1997), we 
calculated for each source the parameter q = ( J — H) — 1.83 x (H — Ks), which determines 
the distance to the reddened vector. YSO candidates, located at the Infrared Excess Sources 
region in Figure 4, have q- values less than —0.15. From this criterium we found 364 YSO 
candidates. The most reddened sources shown in Figure 4 as circled dots present the best 
quality detection in the H and K bands. The numbers correspond to the numbered sources 
in Figure 3. Their locations in the CC diagram indicate that they are likely YSOs, where 
the NIR excess around protostars is due to the optically thick circumstellar disk/envelopes. 
Figure 4 also shows two different groups of sources located between the parallel dashed lines. 
The first one, close to the origin of the diagram, corresponds to several blue sources which 
can be seen in Figure 3. They are probably foreground stars not related to the IRDC, and 
as can be noticed in the CC diagram they are clustered near the main sequence curve with 
visual extinction values from Ay ~ 3 to 8 mag. The second group is located between a 
visual extinction that ranges from Ay ~ 18 to 28 mag. Most of these sources correspond 
to background and possible molecular cloud embedded main sequence stars. 

°| I 1 I 1 I 1 I 1 I 

2MASS 18564827+0118471 - 




H-Ks 

Fig. 4. Color-color diagram of the infrared sources in the vicinity of IRDC 34.776-0.554. The two 
solid curves represent the location of the main sequence (thin line) and the giant stars (thicker 
line) derived from Bessell & Brett (1988). The parallel dashed lines are reddening vectors with the 
crosses placed at intervals corresponding to five magnitudes of visual extinction. We have assumed 
the interstellar reddening law of Rieke & Lebofsky (1985) (Aj /A v =0.282; A H /A V =0.175 and 
Ak/Av=0.112). The plot is classified into three regions: Cool Giants, Normally Reddened Stars 
and Infrared Excess Sources. The most reddened sources are indicated as circled dots and the 
numbers correspond to the numbered sources of Figure 3. 

Figure 5 displays the spatial distribution of the most reddened sources mentioned above 
over the IRAC 8 /im band emission. The emission in this IR band clearly depicts the 
border of the HII region, with the IR radiation mainly originated in the policyclic aromatic 
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Fig. 5. IRAC 8 jj,m band emission obtained from GLIMPSE. The crosses indicate the most red- 
dened sources obtained from the photometric study performed in the area indicated by the box. 
The box represents the same region as shown in Figure 3 

hydrocarbons (PAHs) . The figure also includes the contours of the radio continuum emission 
of the SNR and the HII region and a box that indicates the region where the photometric 
study was performed. From this figure, it is notorious that the most reddened sources lie 
on the border of the HII region that is seen in projection interior to the SNR shell. Also 
the location of these sources delineate the border of IRDC 34.776-0.554. 
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Fig. 6. GLIMPSE color-color diagram, [3.6] - [4.5] versus [5.8] - [8.0], for the most reddened 
sources remarked in Figure 3. Class I and Class II regions are indicated following Allen et al. 
(2004). The ellipse centered at 0,0 encloses the region of main sequence and giant stars. 

Figure 6 shows IRAC color-color diagram for the most reddened sources highlighted 
and numbered in Figure 3. The regions that indicate the stellar evolutionary stages are 
based on the criteria presented by Allen et al. (2004). The sources IR1 and 8 are not 
included because they have no detection in some IRAC bands. Only source 7 lies in the 
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region of protostars with circumstellar envelopes (Class I), the sources 2, 5, 6, 10 and 12 lie 
in the region of main sequence and giant stars (Class III). None source occupies the region 
of young stars with only disk emission (Class II). However the sources 1, 3, 4, 9 and 11, 
located outside the delimited regions, could therefore be reddened Class II objets (Allen 
et al. 2004). 

In particular, in a region of 90" x 90" towards IRAS 18542+0114 (IR1), the region 
analyzed in molecular lines (see next sections), the above 2MASS color criteria (Hanson 
et al. 1997) identifies two sources: 2MASS 18564827+0118471, the most reddened source in 
the whole studied region (box in Figure 5), and 2MASS 18564849+0118597, both indicated 
in the CC diagram (Figure 4). In what follows we analyze the molecular gas in this region. 

3.2. Molecular emission 
3.2.1. HCO+ J=4-3 

Figure 7 shows a 90" x 90" map obtained towards the infrared source IRAS 18542+0114 
in the HCO + J=4-3 line integrated between 40 and 50 km s _1 . As noticed above, in this 
region there are two 2MASS sources, YSO candidates: 2MASS 18564827+0118471 (source 
1 in Figure 7, hereafter IR1) and 2MASS 18564849+0118597 (source 2 in Figure 7, hereafter 
IR2). 

20.0 
10.0 
1:19:00.0 

I 500 
1 

| 40.0 
30.0 
20.0 
18:10.0 

18:56:50.0 49.0 48.0 47.0 46.0 

Right ascension 

Fig. 7. HCO + J=4-3 emission integrated between 40 and 50 km s _1 towards the infrared source 
IRAS 18542+0114. The contours levels are 3.5, 5, 7, 9 and 11 K km s" 1 . The crosses indicate the 
positions of the 2MASS point sources that according to the color criteria described in the text 
could be YSOs. The angular resolution is ~ 22" and the er rms ~ 0.5 K km s _1 . 

From Figure 7 a dense HCO + clump is evident with IR1 exactly lying over the 
HCO + maximum and IR2 towards a border of the clump. By assuming as the limit of 
the HCO + .1=4- 3 structure the 3.5 K km s _1 contour (~ 7 times above the noise level), 
we can estimate a radius of ~ 25" for this clump, 0.36 pc at the adopted distance of 3 kpc. 

Figure 8 shows an HCO + J=4-3 spectrum obtained from the average over the whole 
observed feature. From these observations we derived line parameters of the HCO + J=4- 
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3 that are presented in Table 1. The average main-beam brightness temperature (T m b), 
center line velocity, line width (FWHP), the average integrated HCO + intensity (J T m b dv ) 
and the radius, are given in Cols. 1 to 5, respectively. Errors are formal la value for the 
model of the Gaussian line shape. 

Table 1. Molecular clump parameters from the HCO + J=4-3 emission. 



<T mb ) 


Vlsr 


Av 


(/ T mb dv) 


R 


(K) 


(km s^ 1 ) 


(km s^ 1 ) 


(K km s^ 1 ) 


(pc) 


1.65 ± 0.20 


45.10 ± 0.30 


4.40 ± 0.60 


8.00 ± 1.00 


0.36 



The HCO + J=4-3 profile presented in Figure 8 presents a dip at v ~ 45 km s _1 , the 
central velocity of the whole molecular complex. Only the profiles towards the center of 
the clump, that is towards IR1, clearly display this feature (see Figure 9). Such dip can 
be interpreted as a self-absorption effect caused by less excited gas. As Hiramatsu et al. 
(2007) propose, this kind of spectral features may indicate a significant population of the 
J = 3 level in the outer gas, that requires a moderately high density but with lower J=4-3 
excitation temperature than in the center of the clump. Such effect might be illuminating 
the existence of a density gradient in the clump. 

We have also analyzed the existence of a possible blue asymmetry which would be an 
evidence of infalling gas towards a protostar (Leung & Brown 1977; Zhou 1992). In the star 
formation model of an infalling envelope, the blue-shifted component is stronger than the 
rcdshiftcd one. This is due to the fact that the rcdshiftcd component is absorbed by the 
cooler infalling gas in the near half of the envelope along the line of sight. We did not find 
such an asymmetry in our spectra. However, as discussed by Gregersen et al. (1997, 2000), 
the blue asymmetry line profiles do not appear in the totality of the known YSOs. One 
reason could be that if a rcdshiftcd outflow exists in the telescope beam, the rcdshiftcd 
components absorbed by the infalling envelope would be compensated. 



10 



S. Paron et al.: Star formation activity towards the SNR W44 



20 70 km/s 



if\nAM\K/Ah /It A^lfWir.^-L^.r. 


ftiM i^Vl 1yAl.-4lW1V-1.fi .r 






i i ' | i ' i 


[IrvVVTrVVI 7 | 'If"!/ » 














A 
















1 a At «A \a n . tijA 1 


A 


: A 




■ .. r i n V| . ..h r,,'l .'1 .'■ . 














^vAAaj..,,^, 


■ A 

/^w-^v — 


















A. 



























































20 10 -10 -20 

RA offset ["] 



Fig. 9. HCO + J=4-3 spectra obtained towards IR1, whose offset position is approximately (10", 
0"). The velocity range goes from 20 to 70 km s~ and is shown at the top left profile. Note the 
profiles that present a dip at v ~ 45 km s near the center of the surveyed region. 



3.2.2. 12 CO J=3-2 

Analyzing the 12 CO J=3-2 data we find that the profiles towards IR1 are remarkably 
broader than those towards the border of the surveyed region. For comparison Figure 
10 shows two spectra: the upper one, obtained towards the center of the region and the 
bottom one towards the edge of the map. Such phenomena can be interpreted as evidence 
of an outflow observed along the line of sight driven by the source IR1. In effect, the blue 
component extends from 27 to 37 km s _1 while the red component goes from 50 to 63 km 
s _1 . Figure 11 displays the 12 CO J=3-2 emission integrated in both velocity ranges. The 
blue and red shifted components are presented in thin and thick contours, respectively. 
The crosses indicate the positions of the 2MASS sources IR1 and IR2. 

It is worthwhile to remember that the 12 CO J=3~2 line almost always appears self- 
absorbed towards star-forming regions and often reveals information on the gas kinematics 
(Johnstone et al. 2003). Inspecting our data, most of the 12 CO profiles in the region have 
a dip at v ~ 44 km s" 1 . This velocity is very close to the peak velocity of the 13 CO J=3-2 
line (see subsection 3.2.3), which is an optically thinner line. Such correspondence strongly 
suggests that the dip in the 12 CO profiles is indeed caused by self- absorption by less excited 
gas (see e.g. Zhou et al. 1993). Thus we can discard the possibility of two 12 CO emission 
components with different velocities along the line of sight. 

Figure 12 (left) displays the 12 CO position-velocity diagram along a line at constant 
Right Ascension (18 h 56 m 48 s ), while Figure 12 (right) shows the same diagram along a line 
at constant Declination (+01°18'45"). Redshifted and blueshiftcd features extending in 
velocity are clearly evident in both images. 

Based on these observational results we can confirm that we are observing gas outflows 
originated in the source IR1. These results together with the above noticed slightly extended 
emission in the IRAC 4.5 /im band (indicative of shocked gas by outflows) and its position 
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Velocity [km/s] 

Fig. 10. 12 CO J=3-2 spectra. Upper: spectrum towards the center of the observed region. Bottom: 
spectrum towards the edge of the region. The borders between outflow and line core are shown 
with two vertical lines. The mis noise is a rms ~ 0.15 K 
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Fig. 11. The thin contours represents the CO J=3-2 emission integrated between 27 and 37 km 
s _1 (blueshifted component). The contour levels are 5, 10, 15 and 20 K km s . The thick contours 
show the 12 CO J=3-2 emission integrated between 50 to 63 km s _1 (redshifted component) and 
the levels are 10, 15, 20 and 25 K km s -1 . The crosses indicate the positions of the 2MASS sources 
IR1 and IR2. 



in the infrared color-color diagram (see Section 3.1), confirm that IR1 is a YSO. In Figure 
13 we present a three dimensional picture of the HCO + structure and the outflows mapped 
by the 12 CO. The figure shows the 12 CO and the HCO + data as isosurfaces of intensity. 
The vertical axis corresponds to the velocity while the horizontal axes represent the spatial 
coordinates. The feature at v ~ 40 - 45 km s _1 is the HCO + emission at 0.7 K. The 
extended features in velocity are the 12 CO emission at 1.8 K, which clearly show the 
outflows structure. 
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Fig. 12. 12 CO position-velocity diagrams. Left: along a line at constant RA = 18 h 56 m 48 s . Right: 



along a line at constant dec. = +01°18'45" 





Fig. 13. Three dimensional picture showing isosurfaces of intensity. The feature between ~ 40 
and 45 km s _1 is the HCO + emission at 0.7 K while the extended features in velocity (along the 
vertical axis) show the isosurface of the 12 CO emission at 1.8 K. 

3.2.3. 13 CO J=3-2 and CS J=7-6 

Figure 14 displays the profiles of the 13 CO J=3-2 and CS J=7-6 emission, left and right, 
respectively, obtained towards the center of the studied region (RA = 18 h 56 m 47.8 s , dec. 
= +01°18'45", J2000). The parameters determined from Gaussian fitting of these lines are 
presented in Table 2. T m (, represents the peak brightness temperature, V; sr the central 
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velocity, Av the line width and / the integrated line intensity. Errors are formal la value 
for the model of the Gaussian line shape. 
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Fig. 14. 13 CO J=3-2 profile (left) and CS 3=7-6 profile (right) obtained towards RA = 
18 h 56 m 47.8 s , dec. = +01°18'45" (J2000). The rms noise is cr rms ~ 0.10 and 0.07 K for each 
line, respectively. 



Table 2. Observed parameters of the 1 CO J=3-2 and CS J=7-6 emissions towards RA = 
18 h 56 m 47.8 s , dec. = +01°18'45" (J2000). 



Emission 






Av 


I 




(K) 


(km s _1 ) 


(km s- 1 ) 


(K km s- 1 ) 


13 CO J=3-2 


3.55 ±0.30 


44.10 ±0.20 


5.45 ±0.50 


21 ±2 


CS J=7-6 


0.50 ±0.10 


43.80 ±0.90 


4.00 ±0.40 


2.2 ±0.5 



We have calculated the ratio of the integrated line intensities between the 13 CO J=3-2 
and J=l-0 lines ( 13 R3~2/i-o)- The latter was extracted from the GRS (Jackson et al. 2006) 
and the J=3-2 line was convolved to the J=l-0 beam. We obtain 13 R3_2/i-o ~ 1-4. 

Besides, in order to compare the 13 CO J=3-2 with the 12 CO J=3-2 to obtain an 
isotopic ratio between peak temperatures for the different velocity components, we have 
fit the 13 CO J=3-2 spectrum with three Gaussians. We use the Gaussians with lower and 
higher velocities, closer to the blucshiftcd and redshiftcd components, respectively, as seen 
in the 12 CO emission (see Figure 10). Thus we obtain the following peak temperature 
ratios: 12 T/ 13 T - 13.5 and 12 T/ 13 T ~ 8.8 for the emissions close to the blueshifted and 
redshiftcd components, respectively. Finally, from these ratios we derive the optical depths 
for each line. For this region of the Galaxy, according to Milam ct al. (2005), we use an 
isotopic abundance ratio of [ 12 CO]/[ 13 CO] ~ 70. Assuming an excitation temperature of 
T ex ~ 20 K for both lines, we derive the following optical depths: t 12 ~ 6 and r 13 ~ 0.07 
for the emission close to the blueshifted component and r 12 ~ 10 and r 13 ~ 0.1 for the 
emission close to the redshiftcd one. 
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3.3. Kinematics and dynamics of the outflows 

In this section we investigate the physical parameters of the discovered outflows. To carry 
out these calculations we applied the relations from Choi et al. (1993). The CO column 
density (Np°) and the mass (Mi) for each channel arc calculated from: 

JVJ»= 1.1X10"^^- ^ r and 

D(n,T k ) 1 - exp(-T 32 ) 

where 

D(n,T k ) = f 2 [Ju(T ex ) - J v (T bk )][l - exp(-16.597/T ex )\, 

is the peak temperature of the 12 CO J=3-2 line, Av the channel width, f 2 is the 
fraction of CO molecules in the J=2 state, d = 3 kpc is the distance to the giant molecular 
cloud, Vl is the solid angle subtended by emission. [H]/[C] = 2.5 x 10 3 (Grevesse et al. 1991), 
[C]/[CO] = 8 (Dickman 1978; van Dishoeck et al. 1992), and T32 is the optical depth of the 
12 CO J=3-2 transition calculated above, a value of 1.5 is adopted for D(n,T k ) following 
Choi et al. (1993). With these parameters we have estimated the mass M, the momentum 
P and the kinetic energy E k for each outflow from: 

M = J2 M i: 

i 

P = Y J M\v>-v \, 

i 

Ek=^2^M l (v l -v ) 2 , 

i 

where vq is the systemic velocity, ~ 45 km s" 1 . Finally we calculate the dynamical timcscalc 
tdyn = R/Vch, where V c h = P/M is the characteristic velocity and R is the size of the 
outflows, the mass loss rate M = M/td yn and the average driving force F = P/td yn - As 
the outflows are seen along the line of sight, that is the blue lobe is superposed on the red 
lobe, we estimate their sizes R from the FWHM contour (sec Figure 11). The results are 
presented in Table 3. 

Table 3. Outflows parameters. 



Shift 


R 


M 


tdyn 


M 


P 


E k 


F 




(PC) 


(M @ ) 


(yr) 


(Moyr- 1 ) 


(M km s" 1 ) 


(Melons" 1 ] 2 ) 


(Mokms^yr- 1 ) 


Blue 


0.35 


7.0 xlO -4 


2.6 xlO 4 


2.7 xlO" 8 


7.7 xlO" 3 


4.4 xlO" 2 


3.0 xlO -7 


Red 


0.30 


1.4 xl0~ 4 


4.2 xlO 4 


3.3 xlO" 8 


1.2 xl0~ 2 


5.5 xlO -2 


2.7 xlO" 7 



4. Discussion 

From the 2MASS data study we found that the sources with the largest NIR excess lie on the 
border of the HII region G034.8-0.7 and around of the dark cloud IRDC 34.776-0.554, a sce- 
nario that strongly suggests active star formation. The source 2MASS 18564827+0118471 
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(IR1 in this paper) is the most reddened source in the NIR and is also bright and slightly 
extended in the IRAC 4.5 /im band. This latter emission, which contains H2 and CO vi- 
brational lines, can be produced when protostellar outflows crash into the ambient ISM, 
suggesting that IR1 could be a YSO. 

Our molecular observations have revealed the existence of a HCO + clump where IR1 is 
embedded. It is known that such molecular specie enhances in molecular outflows (Rawlings 
et al. 2004). In effect, a strong enhancement of the HCO + abundance is expected to occur 
in the boundary layer between the outflow jet and the surrounding molecular core. This 
would be due to the liberation and photoproccssing by the shock of the molecular material 
stored in the icy mantles of the dust. Indeed, through the 12 CO J=3-2 line we discovered 
the presence of outflows extended along the line of sight related to the HCO + clump and 
IR1. In this context, we conclude that the source IR1 is a YSO driving bipolar outflows 
which are enhancing the abundance of HCO + . 

The obtained outflows masses are similar to those obtained for the outflows of the source 
IRS6, a low-mass YSO in Cederblad 110 (Hiramatsu et al. 2007). We compared the derived 
parameters for the outflows (see Table 3) with those of the statistical study carried out by 
Wu et al. (2004) based on the properties of a sample of about 391 molecular outflows. In 
that work the sources were divided into outflows with low mass and outflows with high mass 
according to either the available bolomctric luminosity of the central source or the outflow 
mass. Based on the compiled data the authors conclude that the outflow phenomenon is 
common either in low and high mass star formation regions. The red and blue outflows 
masses that we derived for IR1 are slightly smaller than the lowest masses discussed by 
Wu et al. (2004). Therefore, IR1 would correspond to the "low mass group" as catalogued 
by them. It is important to remark that since we made no corrections to the velocities 
for possible projection effects, the momenta and kinetic energies that we derive are lower 
limits. Also, as Choi et al. (1993) note, these effects can increase (if the flow is in the plane 
of the sky) or decrease (if the flow is along the line of sight) the timescale (td yn m Table 
3) by substantial factors. In addition, this time would be smaller if we use the maximum 
velocity rather than V c h- However our tdyn values are very close to the average dynamical 
time of the "low mass group" presented in Wu et al. (2004). Other physical parameters, 
like kinetic energy, momentum and driving force listed in Table 3 also agree with those of 
the "low mass group" . We can therefore conclude that we are probably observing a source 
with outflows of low mass. 

One stage in the formation of a star is the main accretion phase. This is the phase 
during which the central object builds up its mass from a surrounding infalling envelope 
and accretion disk. Observational evidence shows that the main accretion phase is always 
accompanied by a powerful ejection of a small fraction of the accreted material in the 
form of prominent bipolar outflows (e.g. McKee & Ostriker 2007; Bachiller 1996). Class 
protostars are considered the earlier stage of star formation with typical ages of 1 — 3 x 10 4 
years. Most of them, if not all, drive powerful "jet-like" CO molecular outflows, which 
gradually disappear as the protostar goes into the following evolutionary stages. Thus, the 
dynamical times of the outflows can be used as an approximation of the life time of YSOs. 
Therefore we suggest an age of a few xlO 4 years for IR1. 
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Additionally, we performed a Spectral Energy Distribution (SED) fitting of IR1 near 
and mid IR fluxes extracted from the 2MASS Point Source Catalog (J, H and K bands), 
from the Glimpse Catalog (Spitzer-IRAC bands) and from the MSX Point Source Catalog 
(8.28, 12.13, 14.65 and 21.30 fim bands). Based on the SED models 2 presented by Robitaille 
et al. (2006, 2007) and considering a visual extinction between 15 and 30 mag (see Section 
3.1), we suggest that IR1 can be a high mass YSO of ~ 20 M Q . Besides, the models predict 
an age between 7 x 10 4 and 10 5 years for this YSO, and a circumstellar disk of 0.01 M Q with 
an accretion rate of ~ 10~ 6 M Q yr~ 4 . The outflows and accretion rates presented in Table 
3 and derived from the SED models, respectively, are considerably lower than the typical 
rates of massive YSOs (10 -4 M yr _1 in both cases; McKee & Tan 2003). As noticed 
above, IR1 presents a high NIR excess, which suggests that it is embedded in a very dense 
envelope. Taking into account the YSO mass obtained from the SED models, the estimated 
ages and the mentioned outflow and accretion rates for IR1, we suggest that this source is 
probably a massive YSO that could be near the end of the accretion process and has not 
ionized and dissipated its envelope yet. 

The integrated line 13 CO 3-2/1-0 ratio ( 13 R-3-2/i-o) obtained towards the center of 
the region is very similar to that obtained towards a circumstellar disk around a pre-main 
sequence star by van Zadelhoff et al. (2001). As the authors remark, the 13 R.3_2/i-o could 
indicate temperatures of ~ 20 — 40 K. However, care should be taken with the interpretation 
of this result since the emission of the two lines most likely comes from different regions of 
the disk due to the difference in optical depth of the two lines. 

The sulphur-bearing species, like the CS, have been associated with shocks and/or 
outflows. However it is still debatable whether sulphur-bearing species are good diagnostics 
of outflows (Johnstone et al. 2003). According to Moriarty-Schieven et al. (1995), the CS 
J=7-6 line traces high densities (>10 7 cm" 3 ) and warm temperatures (>60 K) in the 
protostellar envelopes. Thus, we propose that the detection of this transition indicates 
the presence of dense and warm gas, showing the inner denser and warmer part of the 
protostellar envelope where IR1 is forming. 

4.1. Could the SNR W44 or the HII region G034.8-0.7 have triggered the formation of IR1? 

From our molecular results and the IR emission study, we conclude that IR1 is indeed a 
YSO. As it lies in the border of the HII region G034.8-0.7, which is in turn evolving in a 
molecular cloud shocked by the SNR W44, we investigate their genetic connection. 

According to the derived outflows dynamical time (see Table 3) and to the SED age 
estimate, IR1 might be a YSO with an age of a few x 10 4 or maybe ~ 10 5 years. Based on 
the age of the associated pulsar, Wolszczan et al. (1991) estimated an age of about 2 x 10 4 
years for the SNR W44. Taking into account that the ages of IR1 and W44 are comparable 
and that a delay between the main triggering agent and the subsequent star formation is 
expected, it is unlikely that the SNR has triggered the formation of IR1. 

On the other hand, as mentioned above, IR1 is seen in the plane of the sky overlapping 
the photodissociation region (PDR) associated with the HII region G034.8-0.7, located at 

2 http:/ /caravan. astro.wisc.edu/protostars 
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the same distance that W44 (Ortega et al. 2007). In the literature there are many works 
about star formation triggered by the expansion of HII regions (e.g. Elmcgrccn & Lada 
1977; Deharveng et al. 2003, 2005; Comeron et al. 2005; Zavagno et al. 2006). In what 
follows we investigate this possibility. From the 12 CO J=l-0 study performed by Seta 
et al. (2004), a dense molecular shell surrounding the HII region G034.8-0.7 is evident at 
the velocity range 35 - 40 km s _1 , which suggests that in this region can take place the 
"collect and collapse" process (Deharveng et al. 2005; Zavagno et al. 2006). In this process, 
the compressed shocked layer generated in the inhomogeneous medium by the expansion of 
the HII region may become gravitationally unstable along its surface on a long timescale. 
The process produces massive fragments of material, which allow the formation of massive 
stars and/or clusters. According to Ortega et al. (2007), taking into account that the 
exciting star of the HII region G034.8-0.7 should be of spectral type between 09.5 and 
03, the age of the HII region would be 10 5 - 3 x 10 6 years. These times are compatible 
with star forming processes on the associated PDR. Besides, as noticed in Section 3.1, the 
most reddened sources in the analyzed region lie preferentially in the border of the HII 
region that is seen in projection interior to the W44 shell. We can therefore conclude that 
an scenario where the expansion of the HII region G034.8-0.7 has triggered the formation 
of IR1 is more possible. Careful modclization of this scenario, including the actions of the 
HII region and the SNR, is planned for the near future. 

Future molecular observations are planned to be performed towards this region in order 
to look for more YSOs embedded in the perturbed molecular cloud. Besides, extended CS 
observations towards IR1 would be important. The transitions of this molecular specie 
are very useful to derive detailed density structure within the molecular envelopes of star- 
forming region. 

5. Summary 

We present an infrared study and molecular observations towards the IR source IRAS 
18542+0114 carried out with the purpose of exploring signatures of star forming activity 
in the border of the HII region G034.8-0.7, which evolves within a molecular cloud shocked 
by the SNR W44. We found that IRAS 18542+0114 presents a slight extended emission 
in the IRAC 4.5 /im band, which suggests that this source may be a YSO driven outflows. 
IRAS 18542+0114 is resolved into several sources in the 2MASS All-Sky Point Source 
Catalog, and two of them, according to color criteria are YSO candidates. One of these 
candidates, 2MASS 18564827+0118471 (IR1 in this work), is the most reddened source in 
the NIR in the studied region. From SED models we suggest that IR1 is probably a massive 
YSO that could be near the end of the accretion process. Our molecular results confirms 
that IR1 is indeed a YSO. These results can be summarized as follows: 

(a) We discovered a HCO + clump towards IR1. It is known that such molecular specie 
enhances in molecular outflows, which arc strong evidences of star forming activity. 

(b) We discovered 12 CO outflows towards IR1 extended along the line of sight. From 
the dynamical time of these outflows and a SED age estimate, we suggest an age of a few 
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xlO 4 or maybe ~ 10 5 years for IR1. The discovery of molecular outflows confirms that IR1 
is a YSO. 

(c) From the detection of the CS J=7-6 line we confirm the presence of high density 
(>10 7 cm~ 3 ) and warm (>60 K) gas towards IR1, probably belonging to the protostellar 
envelope where this YSO is forming. 

(d) By comparing the estimated age for IR1 and the age of the SNR W44 (2 x 10 4 
years), we conclude that it is unlikely that the SNR has triggered the formation of IR1. On 
the other hand, taking into account the age of the HII region G034.8-0.7 and that IR1, as 
other reddened sources, lie on its border, we find that the HII region may have initiated 
the star formation through the "collect and collapse" process. 

Future observations and theoretical modclization are planned in order to further inves- 
tigate the star formation in this complex region. 
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